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Next-to-Leading-Order Staggered Chiral Perturbation Theory * 
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We have extended staggered chiral perturbation theory (SxPT) to Oia^]?), ©(a*), and O(a^m), the orders 
necessary for a full next-to-leading order (NLO) calculation of pseudo-Goldstone boson (PGB) masses and decay 
constants including taste-symmetry violations. We present predictions relating 50(4) taste-breaking splittings in 
masses, pseudoscalar decay constants, and dispersion relations. These can be used to test the v^Det trick. 



1. Introduction 

Lattice simulations with staggered fermionspp 
are fast compared to those with other standard 
fermion discretizations [^15] . They therefore allow 
QCD calculations on a reasonable time scale with 
light dynamical quarks. In addition to this com- 
putational advantage, staggered fermions also 
possess a remnant J7(l) chiral symmetry, even at 
nonzero lattice spacing. Their primary disadvan- 
tage is that each continuum staggered fermion 
flavor comes in four degenerate tastes, and the 
continuum SU (4) symmetry associated with each 
flavor is almost completely broken by the lattice, 
leaving only a discrete subgroup. While taste- 
breaking effects only enter at C'(a^), they turn 
out to be important numerically at current lat- 
tice spacingS(4j- 

In order to extract physical quantities from 
lattice simulations, one must perform both con- 
tinuum and chiral extrapolations using func- 
tional forms determined in chiral perturbation 
theory '5'. Thus, forms applied to staggered sim- 
ulations must account for taste violations. Near 
the continuum limit, where there is a controlled 
expansion in powers of the lattice spacing, SxPT 
systematically describes taste violations in the 
PGB sector. Reference [B] determined the O(a^) 
staggered potential for a single flavor and showed 
that even at tree-level it correctly predicts the 
mass degeneracy pattern observed among the 
different tastes of PGBsTW. References ^ 
and JT] extended the potential to multiple stag- 
gered flavors and calculated the 1-loop mass and 
decay constant for the taste Goldstone (^5) me- 
son. However, it is well known from continuum 



XPT that one must work to at least NLO for 
adequate chiral extrapolation of lattice data|12|. 
Thus, the utility of S^PT is limited without 
higher-order operators. We have therefore deter- 
mined all hnearly-independent O(a^p^), ©(a''), 
and 0(a^m) operators in the staggered chiral La- 
grangian for a general multiple-flavor, partially 
quenched theory, including source terms for left- 
and right-handed currents and scalar and pseu- 
doscalar densities.^ The full NLO staggered chi- 
ral Lagrangian consists of these operators, which 
incorporate discretization effects, plus the stan- 
dard continuum 0{p^^, 0(j?in), and 0{rr?) op- 
erators. 

2. NLO Staggered Chiral Lagrangian 

Determination of the staggered chiral La- 
grangian requires two steps. First one constructs 
the quark level effective action, including ex- 
plicit nonzero a effects. Following the method of 
Svmanzik|14|. S^j f consists of all local continuum 
lattice operators of dimension < 6 invariant under 
the lattice symmetries. In practice it is easiest to 
first determine all allowed lattice operators of di- 
mension < 6, as in Refs. jl5l6| . and then map 
these onto continuum quark- level operators[B]. 
Next one maps the continuum quark-level oper- 
ators onto chiral operators describing the PGB 
sector. To be most general, we consider the par- 
tially quenched (PQ) theory - unqucnchcd QCD 
can be recovered by taking an appropriate limit. 
SxPT uses the following power counting scheme: 
P^I^QCD ~ m/KQCD « a^A^p^, which applies 
to current staggered lattice simulations. Thus the 
LO staggered potential is of 0{a^) 6 10 , while 
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^The full list of NLO taste-breaking operators is given 
elsewhere ll3l . 
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the NLO taste-violating operators in the stag- 
gered chiral Lagrangian are of O(a^p^), 0(a'^), 
and 0{a^m)^. 

The NLO chiral operators can be divided into 
two types according to how badly they break 
the continuum SU (4) taste symmetry. Those of 
type A are rotationally invariant and preserve an 
50(4) subgroup of the continuum SU{A) taste 
symmetry. In contrast, type B operators are only 
invariant under certain combined spin and taste 
rotations, and maximally break the continuum 
symmetry down to the lattice symmetry group. 
Because all operators in the LO chiral Lagrangian 
are of type A, the tree-level PGB masses are split 
into irreps of S'0(4)-taste: ^i,^5,^f^, C^s, and 
However, analytic NLO contributions due solely 
to the B-type operators can cause further mass 
and decay constant splittings among the tastes. 

3. NLO Relations for PGB Properties 

Because there are 234 new NLO taste-breaking 
operators, determining all of their coefficients 
does not seem possible. Nevertheless, there is 
one set of quantities to which only a few con- 
tribute and for which predictions are possible - 
those which are nonzero only because S'0(4)-taste 
symmetry is broken at NLO. 

We first consider S'0(4)-taste and rotational 
symmetry violation in the pion dispersion rela- 
tions. Only eight ©(a^p^) operators contribute^: 

E E {C2Str(9^I]t^^,a^I]e,^) 

+ Cr Str(I]a^Ste^,)Str(St5^Ee.^) 

+ Cio [Str(Sa^E+e^,5]9^St^,^) + p.c] 

+ Ci3 [Str(S9^Et^^^)str(Ea^Et^^^)+p.c.] } 

+ C44A [Str(9^Ste^59^Et^^j 

+ C47V [Str(9^St^^)2-Hp.c.] 

+ C47A [Str(9^St^^5)str(9^Et^^^)+p.c.] }.(1) 

One can see that these operators break 50(4)- 

^Operators of 0{a'^rn) only occur in type A, while type 
B ©(a*) operators contract into S'0(4)-invariant ones for 
two-PGB processes. 



taste and Euclidean rotational symmetry from 
the fact that they have more than two repeated 
indices. This property is common to all B- 
type operators in the staggered chiral Lagrangian. 
Note that four of these operators contain two su- 
pertraces; they produce hairpin correlators which 
only contribute to properties of flavor-singlet 
PGBs. Taste and ^5 hairpins were not present 
in the LO Lagrangian, but do occur at NLO. 

As expected, we observe rotational symmetry 
violation in the pion dispersion relations^: 

Ej = p^l+ Si- Sk) + ml{l+ 64) , 

= ipl + P^) + Pli^+Sk5-S45)+m%{l+Sk5), 

Efm = ipf + pl,)il+6lm-Ski) 

+ Pl + ml^{l+ 5i^), 

Eli - {pf + pL){^+hA-5hn) 

+ pl + ml,{l+5kA). (2) 

Here Ep is the energy of a taste F pion extracted 
from the exponential fall-off of the 2-point func- 
tion along the Euclidean time (4) direction. The 
masses, mp, are the full NLO masses excluding 
the contributions from the type B operators enu- 
merated above, so they fall into irreps of 50(4). 
The corrections, Sp = Srn^/rrijp arise from the 
additional 50(4)-breaking operators, and further 
split the tastes into irreps of the timeslice group, 
the rest frame symmetry group of the staggered 
actionlT^^. As this is the maximum possible 
splitting among tastes on a lattice with zero spa- 
tial momentum, these degeneracy classes hold to 
all orders in SxPT. A testable prediction can 
be made in terms of mass and energy splittings 
among the 50(4) irreps ^^5, and ^^i,: 

El -El ^ ^^pi+p]+^^ 
ml - ml [ml + m4)/2 ' 

ggs - El, ^ ^ ^ p1+P^ + 2pI 

"^fc5 - "^Is (™fe5 + "^l5)/2' 

^Dispersion relations for tastes I and 5 are Euclidean ro- 
tation invariant at this order. 

^Explicit expressions for the Ss are not needed for the 
predictions and are given elsewhere. 1131 . 
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These are now the full NLO masses which are 
directly measurable on the lattice. 

The only Euclidean rotational and 50(4) taste 
symmetry-violating contributions to both PGB 
masses and pseudoscalar decay constants are from 
the same 0{a^p'^) operators, through wavefunc- 
tion renormalization. Thus the splittings within 
50(4) irreps in these two quantities are related: 
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This is the simplest prediction of NLO S^PT that 
one can test on the lattice. It is essential for 
these expressions that the Z-factors are 50(4)- 
invariant, and therefore identical for both tastes 
in the expression. Thus they can be tested using 
bare lattice operators, thereby avoiding ©(a^) am- 
biguities in matching lattice and continuum op- 
erators that could destroy the relationship. Be- 
cause, in the Ss, there are three independent co- 
efficients, and there are three splittings, there are 
no predictions relating mass splittings among dif- 
ferent 50(4) irreps. 

We note that, in the talk at Lattice 2004, we 
presented a different prediction involving the ax- 
ial decay constant: 



( fk-fr \ ^ ( m%-ml \ 



fk+l 



(5) 



and so forth. We now realize that this is incor- 
rect for two reasons. First, it does not properly 
take into account wavefunction renormalization 
when calculating /^. Moreover, even if one cor- 
rectly calculates the contributions from the listed 
0{a^p^) operators, this relationship is broken by 
additional operators containing sources for left- 
and right-handed currents. These operators, in 
which the covariant derivative acts on a taste spu- 
rion, rather than a S-field, affect the axial cur- 
rent, but not the masses or pseudoscalar density. 



4. Conclusion 

We have used NLO SxPT to predict quanti- 
tative relationships between PGB masses, decay 
constants, and dispersion relations. We empha- 
size that, although the underlying lattice sym- 
metries tell us that the approximate degeneracies 
among the 50(4) irreps must be split into the 
true lattice irreps, they do not predict any par- 
ticular relationships among the 50(4) splittings 
in various PGB properties. These come strictly 
from SxPT. These relationships therefore provide 
a simple test of whether the effective field theory 
applies at all to the lattice data, which is based on 
the assumption that the underlying lattice "the- 
ory" using the v^Det trick is not sick. Testing 
them will provide empirical evidence either for or 
against the validity of this trick. 
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